Strong Lensing Constraints on Small-Scale Linear Power 
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We place limits on the linear power spectrum on small scales (k > 50ft Mpc -1 ) using measure- 
ments of substructure in gravitational lens galaxies. We find excellent agreement with the simplest 
ACDM models, and in conjunction with other cosmological probes, place constraints on the neutrino 
mass m„, tilt of the primordial power spectrum n, and mass of the dark matter particle m. We 
find n > 0.94, and for a Harrison- Zeldovich spectrum, m„ < 0.74 eV and m > 5.2 keV, at 95% 
confidence. 



PACS numbers: 14.60.Pq, 95.30.Cq, 95.35.+d, 98.62.Sb, 98.80.-k, 98.80.Cq 
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The amplitude of matter fluctuations on subgalactic 
scales can provide a wealth of information on a wide 
range of physics. For example, the inflaton potential [Q, 
the neutrino mass 0], and the physics of the dark mat- 
ter particles j^, [I| can all cause measurable effects on the 
primordial power spectrum, especially on small scales. 
Hence, measurement of the small scale linear power spec- 
trum is of great interest. Unfortunately, direct probes 
of the linear power on these scales has proven difficult, 
due to a host of astrophysical processes which compli- 
cate the interpretation of attempted measurements. The 
most promising avenue of attack has been study of the 
Lyman a forest absorption power spectrum || , however 
it has been argued that inference of the linear power 
spectrum from the Ly a forest will be difficult due to 
nonlinear effects ||. One relatively clean probe of the 
linear power spectrum is the number of collapsed halos 
as a function of mass. Assuming Gaussian density fluc- 
tuations, simple smoothing arguments may be applied 
to calculate the collapsed halo abundance given the lin- 
ear power spectrumQ. This Press-Schechter formalism 
reproduces the mass function produced in cosmological 
simulations remarkably well, with only minor modifica- 
tions [pi] . A consequence of Press-Schechter theory, subse- 
quently confirmed by numerical simulations, is that sim- 
ple inflationary cold dark matter (CDM) models predict 
halos with significant amounts of satellite substructure 
||. The dearth of observed Local Group dwarf satel- 
lites, relative to these predictions, could in principle hint 
at suppression of fluctuations on small scales, implying 
new physics in inflation [jj] or in the dark matter sec- 
tor H U . Unfortunately, once again uncertainties in the 
astrophysics of such dwarfs, such as star formation or 
feedback processes 0, limit the utility of these objects 
in constraining the matter power spectrum. Instead, a 
direct probe of the satellite mass (as opposed to the light) 
is required. Gravitational lensing, as an effect sensitive 
to mass rather than light, can provide such a probe. Re- 
cently, the abundance of satellite substructure in a sam- 
ple of elliptical galaxies has been measured using gravi- 
tational lensing jy]. In this Letter, we quantify the con- 
straints of this detection on the linear power spectrum. 



Gravitational lens galaxies require a significant frac- 
tion, 0.006 < / sat < 0.07, of their mass to be in subhalos 
in order to be consistent with observations of lens im- 
age fluxes and relative positions JO). The most likely 
candidates for such satellites are the CDM subhalos pre- 
dicted by Press-Schechter arguments and seen in N-body 
simulations; see Ref. 11 for more details. Using the 

121 , which has been shown to 



conditional mass function 
match well the satellite mass function found in high reso- 
lution N-body simulations jj] , this / sat directly translates 
into a measurement of the linear rms density fluctuations 
on such mass scales, <r(M sa t). For a Sheth-Tormen (ST) 
multiplicity function flq], we have 



/ sat (> M sat ) = A 
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Here, 8 C — 1.68 is the critical linear overdensity for 
collapse flijfl , erfc(x) is the complementary error func- 
tion, T(n,x) is the incomplete gamma function, A^ 1 = 
l + 2~ p 7r~ 1 / 2 r(i— p) is a normalization constant, D(z) is 
the linear growth factor normalized so that D{z = 0) = 1, 
and Zg a i « 1 is a typical redshift at which the lens galaxy 
halo of mass M ga i ss 10 12 ' 5 M Q formed. For a Press- 
Schechter multiplicity function, we would have a = 1, 
p = and A = i, however ST find that N-body simula- 
tions are better fit by a — 0.707, p — 0.3 and A w 0.322 
M. We estimate the typical formation redshift for satel- 
lites by assuming they are tidally truncated and set- 
ting their inferred density equal to 200 times the back- 
ground matter density when they formed. Assuming a 
flat ACDM cosmology with fijf = 0.33 and Hubble con- 
stant h = 7Jo/(100 km/s/Mpc)=0.66, and adopting the 
ST fit, the lens constraints on satellite Einstein radius 
and mass fraction translate into constraints on a(M). 
This is plotted in Figure [l]. 

Having determined the rms linear mass fluctuations on 
subgalactic scales, we now compare this result to the pre- 
dictions of various cosmologies. A representative sample 
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fi^ft. =0.005- 
0^=0.015- 
m= 10 keV - 
m=3 keV 
n=0.9 



since 



■m v = 9ifl„h 2 eV. 
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FIG. 1: Variance a on mass scale M. Contours depict the 
68%, 90% and 95% confidence regions from lensing. The 
curves depict predictions for the cosmologies denoted which, 
unless labeled otherwise, are flat ACDM universes with h = 
0.66, n M h 2 = 0.15, Q v = and n = 1. For the WDM models, 
m denots the WDM mass. 



are plotted in the figure. Fortuitously, the wide degener- 
acy in the lens constraints happens to lie nearly orthog- 
onally to the direction in which different <r(M) curves 
shift as the cosmology is varied. Clearly, the lensing 
measurement can strongly constrain a variety of physics. 
We consider three such constraints: the tilt of the pri- 
mordial power spectrum, the mass of the neutrino, and 
the mass of the dark matter particle. Throughout, we 
use COBE normalized power spectra fl4|| , assume a flat 
ACDM cosmology with tt DM h 2 = 0.13, tt B h 2 = 0.02 and 
h = 0.66 ± 0.05, and employ fitting functions for the 
transfer functions || fL5f . 

The tilt of the primordial power spectrum, n, is related 
to the slow roll parameters of the inflaton potential [jl6| . 
Recent measurements |l7[] of the cosmic microwave back- 
ground radiation constrain n w 1± 0.1 (1 a) on large 
scales ranging up to the horizon size. Lensing probes 
mass scales some 18 orders of magnitude smaller, so it 
provides the strongest constraints on the tilt that are 
currently possible. We find 0.93 < n < 1.04 at 95.5% 
confidence, consistent with a scale-invariant Harrison- 
Zeldovich spectrum of perturbations. If we instead de- 
rive one-sided limits, we find n > 0.94 at 95.5% confi- 
dence. Inflationary models with broken scale invariance 
Jl, 18 arc tightly constrained; the inflaton potential can- 
not have features which break scale invariance on scales 
larger than ~ 20 kpc. 

Experiments probing atmospheric and solar neu- 
trinos |2(| have established that neutrinos have mass 
and thus constitute particle (hot) dark matter. How- 
ever, these experiments do not provide a definite mea- 
sure of the mass of neutrinos, but instead the squared 
mass difference A(m 2 ) between different neutrino flavors. 
The linear power spectrum constrains the neutrino mass, 



Note that this relation only applies for active neutrino 
flavors, and for vanishing lepton asymmetry, however re- 
cent work indicates that the asymmetry (and hence its 
corrections) are small [^l] . Increasing the energy density 
in neutrinos has the effect of suppressing power on small 
scales. Previous constraints Q from Lyman a forest ob- 
servations limit m v < 5.5 eV. By probing even smaller 
scales, lensing further tightens this bound. We explored 
a range of neutrino masses ranging from 0.05 eV, near the 
lower limit imposed by atmospheric neutrino oscillations 
|l9| up to 3.2 eV. For a Harrison- Zeldovich spectrum, and 
assuming a uniform prior on m u , we find m v < 0.74 eV 
at 95.5% confidence. If we instead assume a logarithmic 
prior, which may be more appropriate, we find m u < 0.53 
eV at 95.5% confidence. Interestingly, the lensing upper 
limit on the sum of neutrino masses is quite close to the 
lower limit on neutrino mass measured by the LSND ex- 
periment |22j, roughly V Am 2 ~ 0.3 eV. 

Additionally, our measurement provides a lower limit 
on the energy spectrum of dark matter particles. We 
follow standard convention and characterize the energy 
spectrum by the mass of a neutrino with an equivalent 
free streaming scale. Warm dark matter models, with 
m < 1 keV, have been proposed as a modification to 
CDM, in order to account for the dearth of dwarf satellite 
galaxies and the lack of cusps in dark matter density in 
the central regions of galaxies The close agreement 
of lensing observations with CDM predictions indicates 
that the free streaming scale must be below the ~ 20 kpc 
scales probed by lensing. Assuming a uniform prior in the 
range 1 keV < m < 32 keV, and assuming a Harrison- 
Zeldovich spectrum, we find that m > 5.2 keV is required 
at 95.5% confidence. Since the dark matter particle mass 
can be much larger than the range we have explored, our 
bound is conservative. 

Certain effects which we have not considered will tend 
to suppress the effects of substructure, and hence bias our 
measurements to low a. First, the conditional mass func- 
tion ignores the effects of tidal stripping and disruption 
of satellite subhalos. Secondly, the finite angular sizes of 
radio QSO's tend to wash out the lens effects of halo sub- 
structure. A proper treatment of these effects is beyond 
the scope of this Letter, however we note that both ef- 
fects tend to increase the required abundance of satellite 
substructure and strengthen the constraints. Hence our 
limits should be regarded as conservative. 

These constraints should all improve as the sample of 
suitable lens systems expands, and as the quality of the 
data improves. In general, multiply- imaged compact ra- 
dio sources will be the preferred probe because of the 
high imaging resolutions available for radio sources using 
Very Long Baseline Interferometry (VLBI). For the cur- 
rent lenses, deep VLBI observations of extended, thin jet 
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FIG. 2: Likelihood of (a) tilt of primordial power spectrum, (b) neutrino mass, (c) dark matter particle mass. 



structures in lenses |23fl are a promising means of break- 
ing the wide degeneracy between satellite mass and mass 
fraction. The Extended Very Large Array (24|] and the 
proposed Square Kilometer Array radio array p5[ can 
easily expand the sample of radio lenses by more than 
an order of magnitude, with a corresponding reduction 
in the uncertainties. 
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